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Abstract 

The predictions of the theories of solid-particle erosion of brittle 
materials are compared to experimental results of studies in which angular 
AI 2 O 3 particles with mean diameters D of 23-270 un are used to erode (111) 
surfaces of silicon single crystals at impact angles a from 20-90° and 
velocities v from 30-150 m/s. The description of the steady-state erosion 
rate by a power law, /W « (v sina)°D™ must be modified to include threshold 
and plasticity effects. Rirthermore the velocity exponent n depends on 
D. Results using abrasives of different sizes mixed together can be 
explained using a logarithmic-normal distribution. The results of 
transient experiments can be used to explain the synergistic effects which 
are observed using a biomodal distribution of abrasives. 


I. Introduction 


The erosion of materials by solid-particle impacts is an important 
process which may limit the service lifetime of components. Brittle 
materials have potential uses in many high-technology energy applications, 
e.g. valves in coal gasification plants, gas turbine blades, electrodes and 
regenerative heat exchangers for MHD applications, and photovoltaic 
devices. Therefore, understanding the erosion process in brittle materials 
is important. This paper will review the progress made in the last two 
years in understanding the erosion process in silicon single crystals, a 
material which not only has applications as photovoltaic devices, but 
represents an ideal brittle solid and, therefore, is important as a model 
material that should closely conform to theoretical predictions. 


II. Theory 

Material removal by impacting particles occurs by lateral crack 
formation i.e., subsurface cracks parallel to the impacted surface, which 
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form primarily as a result of residual elastic-plastic stresses under a 
sharp indentor. Two models based on this experimental observation have 
been proposed to describe the erosion process in brittle materials. Both 
theories assume that the material removed is given by the area containing 
the lateral cracks times the depth of the lateral cracks, which is assumed 
to be proportional to the depth of penetration of the impacting particle. 
Both models assume that the lateral crack size c is proportional to the 
size of the radial cracks, l.e., cracks normal to the impacted surface, 
which form as a rosult of elastic-plastic loading stresses under a sharp 
indentor. In turn, the latter may be viewed as end-loaded half-penny 
cracks, the loading being due to the plastic zone expansion, where the 
plastic zone size is small compared to.the final crack arrest size. 

Fracture mechanics gives c “ (P /K for this situation with P_„„ 

being the maximum contact force and K^. the fracture toughness. The 
theories differ, however, in the calculation of contact stress 
Po “ Pmax/^ax* 4ax contact area. 

The quasi-static model of Wiederhorn and Lawn^^^ calculates the force 
based on the conversion of the kinetic energy of the impacting particle 
modelled as a sharp indentor into plastic work. On the other hand, the 
model of Evans, et al.^^^ neglects plasticity and the contact pressure is 
assumed equal to the dynamic pressure when a spherical particle first hits 
the surface. The depth of penetration is determined from the time of 
contact, and the mean interface velocity, both of which are calculated from 
a one-dimer Honal impact analogue. Both models predict that the steady- 
state erosion rate (weight loss [gl/total weight of abrasive impacting [g]) 
is given by « R™v'^, where R is the particle radius and v is the 
velocity. The exponents predicted on the basis of the two models are given 
in table 1. It may be seen that the only way to distinguish between the 
models lies in a determination of the velocity exponent, n. 


TABIJE 1. Predictions of erosion models. 


Model 

Particle 

Shape 

m 


n 


Quasi- Stat ic^ 

Sphere 

Angular 

2/3 

2/3 

(.67) 

(.67) 

11/6 

22/9 

(1.8) 

(2.4) 

Pulse-Impact 

Sphere 

Angular 

2/3 

2/3 

(.67) 

(.67) 

19/6 

5 

(3.2) 


III. Experimental 

Angular AI 2 O 3 particles were used to erode (111) Si single crystals 
using a slinger-type device.'^' The experimental details have been 
described previously. Single Impacts are examined using scanning 
electron microscopy (SEM). Erosion rates at a fixed impac*" angle. 
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velocity, and particle size are determined from sequential weight-loss 
measurements. 


IV. Results and Discussion 
IV- 1 . Single Impacts 

A typical SEM of a single Impact produced at a = 90° and v = 108 m/s 
using 270- un AI 2 O 3 is shown in Fig. 1. Each fan that orignates from the 
impact site is formed by propagating lateral cracks which periodically 
diverge up to the free surface, causing material removal. The lateral 
crack formation is considered in detail by Evans, et al.^'^^ High 
dislocation densities under the impact sites have been observed,'^’' and 
this provides .^vide’ice for the importance of plasticity. 

Rirther impacts produce overlapping damage sites until eventually a 
steady-state AW is achieved. Figure 2 illustrates a weight-loss curve 
measured for a= 90°, v = 108 m/s, using 37-wn particles to erode a surface 
previously eroded into steady state using large 270- wn particles. The 
erosion rate (the slope) initially decelerates as opposed to an 
accelerating ZW which is always observed on pristine surfaces. The shape 
of the transient is therefore determined by the initial condition of the 
surface. 



Fig. 1. (Left) SEM of single impact produced using 270- un AI 2 O 3 at 
V = 108 m/s and a ,= 90°. 

Fig. 2. (Right) Weight loss as a function of dose for 37- im particles 
impacting a surface previously eroded into steady state (using 
270-iJn particles) at v » 108 m/s and a = 90°. 

IV-2. Particle-size Dependence 


The particle-size exponent m is close to the 2/3 predicted by the 
models for large particles. However, the relation greatly overpredicts ZW 
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for small particle sizes. This indicates that the expression must be 
modified to allow for threshold effects which seem to be manifest for 
smaller particles. If the data for erosion rate ^ and particle size R of 
ref. 4 is plotted as (to Af)/( l-R/R^)^ vs InCR-R^)^'^ to allow a comparison 
to the models that predict the volume removed per number of Impacts, it is 
found that a velocity-dependent threshold size ^ can be obtained such that 
the relationship vR^^'^ •• (1280 ± MO) x 10“® mVs is approximately 
valid. Fbr a* 90^, v >• 100 m/s, the threshold R^ * 6 un. The threshold 
can be related to a critical force required to propagate a crack, and the 
quasi-static model' predicts vR^^'^ ■ constant, while the pulse-impact 
model predicts vR^ ” constant. While the exact relation is difficult to 
evaluate, it appears from the data that the former relation is more 
reasonable. 

IV-3. Particle-Size Distribution Effects 

One of the difficulties %rhich arise in the appraisal of threshold 
effects is that OH depends on the particle-size distribution.^”' Figure 3 
shows the effect of particle-size distribution at a ■■ 90° and v - 100 m/s, 

VELOOTY Im/sl 




InUlmAl) 

Fig. 3. (Left) The normalized erosion rate (measured steady-state 
rate/rate for o - o) as a function of the particle-size 
distribution a. 

Fig. 4. (Right) The logarithm of the steady-state erosion rate as a 

function of the logarithm of velocity for particle sizes of 23, 
37, 130, and 270- im. 
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as measured by the width of the distribution o, on the normalized erosion 
rate (measured rate/rate for o). The solid curve that is calculated 
from a logarithmic-normal distribution is seen to adequately describe the 
experimental results. 

IV-A. Vel o '- 1 - d ependence 

The depend;' t'e of l in v has been systematically Investigated. 

The logarithm of the steau/-state erosion is plotted as a function of the 
logarithm of velocity for four particle sizes in Fig. 4. Good fits to 
[U « v" are obtained, but there is a dependence of the velocity exponent on 
the particle diameter D as shown in Fig. 3. The figure also shows data 
obtained on two different types of silicon carblde'®*^^ which show that Si 
is not unique li> this respect. 

The velocity exponent varies from 3.4 for 37-un particles to 2.55 for 
270-un particles. As seen from Table 1, no current model can explain this 
variation of n with D. It may be postulated that smaller particles have 
shorter contact times and therefore must be approximated using the pulse- 
impact model, while larger particles more nearly satisfy the quasi-static 
model. This predicts a trend in the direction observed. It is interesting 
to note that the velocity exponent obtained using large (1.58-mm diameter) 
spheres impacting MgO is *> 2.1,'^®' in agreement with the trend predicted 
in Table 1. It is believed thaw the exponent for hot-pressed SiC is low 
because of the presence of weakened grain boundaries'^^ which affect the 
erosion Tate. In this respect polycrystalline MgO behaves like MgO single 
crystals' probably because the polycrystalline MgO was relatively pure. 

IV-5. Angular Dependence 

The power-law expression for tU is valid only for normal incidence for 
which is maximum for a brittle solid. For oblique impact angles the 
velocity V can be resolved into a normal component v slnaand a tangential 
component v coscu If frictionless contact conditions exist, only the 
normal component contributes to the erosion, and it can then be given by 
4M (v sino)*'. Normalized data ( iW( o)/ ^(90®)) obtained for various 
velocities and particle sizes are shown in Fig. 6, where the solid line 
denotes 8in^’®ou 


'..te assumption that the tangential component of v does not contribute 
to Hi breaks down for a < 45® where the actual losses are 2-4 times greater 
than those predicted by the model. The additional contribution to /W for 
smaller a can be rationalized if it is assumed to be due to the tangential 
velocity component that arises because of a plastic-deformation cutting 
process, which in a ductile material has a maximum for a •> 20®.'^^^ This 
is consistent with TEM observations'^^ which Indicates that plasticity 
contributes to the erosion process. 


IV-6. Synergistic Effects 

The experimental conditions used in these studies cover the range of 
particle sizes, velocities, and impacts generally expected in service 
applications tihere the components are subjected to an erosive environment 
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l.e. photovoltaic devices unprotected from a dust environment. The models, 
however, can not be assumed adhoc to apply to complex service conditions 
where, for example, several particle sizes or velocities are present 
simultaneously. The simplist assumption Is to use a principle of linear 
superposition which requires that the damage processes occur Independently 
of each other. This assumption Is not In fact valid, and has recently been 
examined' in detail. 

Figure 7 presents the results of an experiment designed to examine 
linear superposition for erosion using a mixture of two sizes of 
particles. The steady-state erosion rate In Fig. 7 Is plotted as a 
function of the weight fraction of the 270- (Si particles (f 27 n) ^ mixture 

of 37- un and 270- wn particles. The simple "law of mixing" given hy tU * 
^270^270 ^37^37* ''here foy is the weight-fraction of 37-un particles 

and the tU's are the respective steady-state erosion rates obtained for 
that size of particles, Is shown as the dashed line and Is not a valid 
description. 


0[^m] 


23 37 (30 270 



ln(D[^m]) 


Fig. 5. (Left) The logarithm of the velocity exponent n as a function 
of the logarithm of the mean particle diameter D. 

Fig. 6. (Right) The steady-state erosion rate normalized by for 
a- 90° as a function of Impact angle a. 

As can be seen from Fig. 2 the erosion rate of the 37-un particles 
impacting a surface pre-eroded with 270- un particles Is initially enhanced 
over the eventual steady-state rate. The enhanced (initial) erosion rate 
tH' 27 he used' describe the results using tU • ^270^70 
^37 *37 + ^270^^ ~ ^270^^ ^'37 ~ ^7^ where the superscript o Is used to 
denote the rate for particles acting Individually and the subscript denotes 
the size. This relation, which requires an accurate measurement of the 
transient erosion rate ^'371 is shorn by the solid lines In Fig. 7 for 
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6 V *> 8 , the range estimated experimentally. The results 

support the predicted trend. It should be mentioned that In an actual 
service application the situation Is most likely to be more complex, due to 
more complicated particle distributions. 



'm 


Fig. 7. The steady-state erosion rate obtained for mixtures of 37 ano 
270-vra particles shown as a function of weight-fraction of 
270- un particles. The dashed and solid lines are explained In 
the text. 

V. summary 

The existing models adequately predict the functional dependence, on 
velocity and size of Impacting particles, of the steady-state erosion rates 
In SI single crystals measured using angular AI 2 O 3 p.rrtlcles If they are 
modified to Include: 1. plasticity for small Impact angles, 2. 

partlcle-slze (and possibly velocity) threshold effects, 3. a partlcle- 
t .ze dependent velocity exponent, and 4. a partlcle-slze distribution 
effect. The above effects are known to exist, but further systematic 
experiments are needed to establish the phenomenology in other systiJis, and 
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to provide a sound basis for the proper relationships needed In physical 
models. Theoretical work Is needed to Incorporate these effects Into the 
models. Synergistic effects are known to exist, but our understanding of 
them Is not complete, and It Is certainly not possible to predict complex 
synergistic effects on the basis of our current knowledge. Finally, the 
projecticle properties (shape and hardness) have never been Investigated. 
Hicrostructural effects in polycrystalline SI are also possible. 
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DISCUSSION: 


BOUJIKIAN: On your formula where DO was 12 microns — at Chat point your erosion 

sort of stopped — are the particles interfering with each ocher? 

ROUTBORT: No. We are extremely careful to feed very slowly so Chat we get no 

interference of particles, in fact single streams of particles. What %re 
cannot do of course, is to get particles 10 microns in diameter to erode. 
What you can do, however, is to calculate what the theoretical threshold 
should be and it turns out Chat in the case of silicon, the threshold is 
even less Chan we can measure in velocity. We have Co go down to 10 
meters per second. There the erosion rate is so slow that we couldn't 
measure it. So it's not a particle interference effect, it's probably a 
real threshold effect, but we haven't proved it unambiguously. The 
material removal rate is proportional to the particle size Co the 2/3 
power, and Che velocity anywhere from Che 2nd to Che 4th potrer. It 
depends on particle shape, because that depends on Che contact condi- 
tions. It depends on the hardness of the material. It depends on what- 
ever model you use, it depends on the acoustic impedance of Che target 
compared to the particle and it depends on the density of Che impacting 
particle. 

BOUJIKIAN: The hardness of the particle, therefore, comes into it. 

ROUTBORT: The hardness of Che target, not Che particle. Yes. 

CHEN: Your model is based on Che complete brittle fracture model, brittle 

material, no plasticity occurred during Che impact. 

ROUTBORT: No. That's not quite right. Because you do assume that Che kinetic 

energy of the indentor, if you will, is converted to plastic work in the 
plastic zone. 

CHEN: No, I'm referring to the Weiderhorn paper about six months ago. He 

used a high-speed camera, and shooting Che particle on Che glass surface, 
he definitely showed there's a scooping. Showed the particle really 
pushed into Che glass surface, and melted it... with energy so high it 
melted the surface and scooped part of the material out.RO 

ROUTBORT: Many people observe intense shear zones where there's actually molten 

material. We have never observed it in silicon. You can indeed calculate 
that there's enough kinetic energy of the impacting particle to melt the 
material depending on the conductivity of Che material. But we've never 
seen it. 

HEIT: Are the abrasive particles directed against the work in an airstream? 

ROUTBORT: No. That's not an airstream. It's under vacuum. ..Che whole system. 

It's under vacuum because this arm is rotating at 10,000 rpm, and it 
doesn't rotate very well in air. The particles are mechanically 
accelerated out Che end of the Cube. 

HEIT: How do you determine the weight loss of the silicon? 



ROUTBORT: We do it sequentially, we put in a charge of 10 grams, tre erode 

away, we stop, we open the vacuum, we take the samples out and weigh them, 
we put them back in. 

HEIT: Is there any embedment of the abrasive? 

ROUTBORT: Absolutely none. We've used dispersive X-ray analysis and there's 

no trace of aluminum; silicon yes, but none of aluminum. Alumina we find 
embedding in all of our metal «K>rk. In fact, many of our metal samples 
that we've run for various reasons or other gain weight due to embedding. 
The aluminum and nickel are both fairly soft, the abrasive particle is 
very sharp, it just sticks in. 

WOLFE: I want to congratulate you on a marvelous piece of work that really 

helps to illuminate what is going on in this silicon removal area. You 
recommended looking at a higher-density particle like aluminum oxide 
rather than silicon carbide or such. I think there's probably a small dif- 
ference. What you did is probably directly applicable to the something 
like sandblasting, while what ve have is a backup of the particles with 
the tool, so the tool actually imparts the velocity onto the particle and 
so therefore the density of the particle is probably not as important as 
its hardness. I think the hardness comes in the size of the impact area. 

If you have a more ductile particle impacting, the impact area is probably 
larger, because the particle spreads out. When you have a very hard parti- 
cle, the impact area is smaller, we have a larger force on a smaller area. 

I suspect that goes more rapidly to the cutting rate question than the den- 
sity in this type of cutting we are doing here. 

ILES: This is the first paper we've had where people are discussing the me- 
chanics of erosion. It seems to me we've got liquid drops and we've also 
got particles of silicon from the kerf, coming at very high speed, loose, 
not bound on the diamond wheel. Are we in the range of speeds where we 
would expect to see some impact with silicon by silicon itself, which would 
perhaps modify the cut rates? 

ROUTBORT: Do you have any idea what the velocities are? 

ILES: I suspect it's in the range of 100 meters per second. 

ROUTBORT: We have significant losses at 10 meters per second with hard 

particles. 

ILES: I'm glad your talk opened up that sort of possibility. That's very 

interesting. 

WOLF: Danyluk's experiments seem to indicate, in light of what you have been 

showing us here now, that depending upon what kind of lubricant ve are 
using, we could get predominantly ductile erosion, or predominantly 
Lrittle erosion. Possibly one kind of hammering of the particles due to 
some tool vibration and so on, and the other kind, just pushes ductily the 
SMterial away. Maybe we can learn to take advantage of these. 

ROUTBORT: These things make a difference of a factor of 4 or so in erosion 
rate. At least the stuff we've studied. Now four is evidently enough for 
you people to make big savings. 



